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Abstract 

(N 

' In the present paper, the phenomenological model developed by the authors in 

OO ■ previous papers has been used to evaluate the degradation induced by hadron ir- 

\ radiation at the future accelerator facilities or by cosmic protons in high resistivity 

^ \ silicon detectors. The damage has been analysed at the microscopic (defects produc- 

ed ' tion and their evolution toward equilibrium) and at the macroscopic level (changes 

, in the leakage current of the p-n junction). The rates of production of primary de- 

i-C \ fects, as well as their evolution toward equilibrium have been evaluated considering 

explicitly the type of the projectile particle and its energy. Vacancy-interstitial an- 
Oh. nihilation, interstitial migration to sink, complex (VP, VO, V2O, CiOi and CiCg) 

(— I \ and divacancy formation are taken into account for different initial silicon material. 

The influence of these defects on the leakage detector current has been calculated 
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1 Introduction 



The use of silicon detectors at the new generation of accelerators or in space 
experiments poses severe problems due to changes in the properties of the ma- 
terial after long time irradiation, and consequently influences the performances 
of detectors. 

The phenomenological model developed by the authors in previous papers 
has been used to evaluate the degradation induced by hadron irradiation at 
the future LHC accelerator facilities and by cosmic protons in high resistivity 
silicon detectors, in two types of silicon: "standard" material (10^^ cm~^ atoms 
of phosphorus, 2 ■ 10^^ cm~^ atoms of oxygen, and 5 ■ 10^^ cm~^ atoms of 
carbon), and "oxygened" one (containing 10^^ cm~'^ atoms of phosphorus, 
4 ■ 10^^ cm~^ atoms of oxygen, and 5 • 10^^ cm~^ atoms of carbon) respectively. 
The damage has been analysed at the microscopic (defect production and 
their evolution toward equilibrium) and at the macroscopic level (changes in 
the leakage current of the p-n junction). 

These theoretical estimates permit to draw conclusions about the damages 
induced in silicon in various irradiation conditions and for different semicon- 
ductor materials. 



2 Main hypothesis of the model 

In the model, the effects of irradiation conditions and various initial impurities 
in the starting material are discussed in the a quantitative manner that the 
defect production and their evolution toward stable defects during and after 
irradiation in silicon is calculated. The model supposes three steps. 

In the first step, the incident particle, having kinetic energy in the intermedi- 
ate up to high energy range, interacts with the semiconductor material. The 
peculiarities of the interaction mechanisms are explicitly considered for each 
kinetic energy [1,2]. 

In the second step, the recoil nuclei resulting from these interactions lose 
their energy in the lattice. Their energy partition between displacements and 
ionisation is considered in accord with the Lindhard theory (see reference [3] 
and authors' contributions [4]). 

A point defect in a crystal is an entity that causes an interruption in the lattice 
periodicity. In this paper, the terminology and definitions in agreement with 
M. Lannoo and J. Bourgoin [5] are used in relation to defects. 



2 



We denote the displacement defects, vacancies and interstitials, as primary 
point defects, prior to any further rearrangement. After this step the concen- 
tration of primary defects is calculated. 

The concentration of the primary radiation induced defects per unit fiuence 
[CPD) in silicon has been calculated as the sum of the concentration of defects 
resulting from all interaction processes, and all characteristic mechanisms cor- 
responding to each interaction process, using the explicit formula (see details, 
e.g. in references [6,7]): 

CPD {E) = ^ fy:(^] L{Eni)s.dn = —^^^^NIEL (E) (1) 

where E is the kinetic energy of the incident particle, Nsi is the atomic density 
in silicon, Asi is the silicon atomic number, Esi - the average threshold energy 
for displacements in the semiconductor, Em - the recoil energy of the residual 
nucleus produced in interaction z, L^Eni) - the Lindhard factor that describes 
the partition of the recoil energy between ionisation and displacements and 
{da/dQ)i - the differential cross section of the interaction between the incident 
particle and the nucleus of the lattice for the process or mechanism i, respon- 
sible in defect production. Na is Avogadro's number. The formula gives also 
the relation with the non ionising energy loss (NIEL), the rate of energy loss 
by displacement {dE/dx)ni [8,9]. 

The kinetic energy dependence of CPD for pions (from reference [10]) and 
for protons [8,11] is used in the present calculations. CPD versus the kinetic 
energy of pions and protons is presented, e.g., in Figure 1 of reference [12]. 

The basic assumption of the model is that primary defects, vacancies and 
interstitials, are produced in equal quantities and are uniformly distributed in 
the material bulk. They are produced by the incoming particle, or thermally 
- only Frenkel pairs are considered. The generation rate of primary defects is 
a sum of two components: 

G = Gr + Gt (2) 
where Gr accounts for the generation by irradiation, and is calculated as: 

Gr = GPD{E) X (^{E) (3) 

with $(i?) the flux of considered incident particles, and Gt for thermal gen- 
eration. 

In silicon, vacancies and interstitials are essentially unstable and interact via 
migration, recombination and annihilation or produce other defects. The con- 
centration of primary defects represents the starting point for the following 
step of the model, the consideration of the annealing processes, treated in the 
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frame of the chemical rate theory. A review of previous works about the prob- 
lem of the annealing of radiation induced defects in silicon can be found, e.g. 
in Reference [13]. 

After silicon irradiation, the following stable defects have been identified ex- 
perimentally and confirmed: Si, VP, VO, V2, V2O, CiOi, Ci, CiCg (after the 
compilations from references [5,14]). Vacancy-interstitial annihilation, intersti- 
tial migration to sinks, divacancy and vacancy - impurity complex formation 
{VP, VO, V2O, Ci, CiOi and CiCg) have been considered supposing the fol- 
lowing chemical reaction scheme: 

V -\- 1 ^annihilation (4) 

G 



I ^ sinks (5) 



VO is the A centre. 



VP is the E centre. 



y + O g (6) 



V + P^VP (7) 



K3 

V + V ^V2 



V + VO^ V2O (9) 



/ + C,^Q (10) 



A + l'^O 



:iii 



(12) 
(13) 



Ki with i = 1, 3-^8 are the reaction constants. Some considerations about the 
determination of the reaction constants are given in references [12,13] and their 
concrete values are given in reference [15], and [16]. Without free parameters, 
this model is able to predict the absolute values of the concentrations of defects 
and their time evolution toward stable defects, starting from the primary 
incident particle characterised by type and kinetic energy. 



4 



3 Radiation environment at the LHC accelerator and in the near 
Earth orbits 

In the present paper, two types of appHcations of silicon detectors are empha- 
sised: for the tracker of experiments at the LHC accelerator, and for space 
missions in the near Earth orbit, as, for example experiments at the Interna- 
tional Space Station. 

At the luminosity of 10^^ cm~^s~^, and assuming an inelastic non-diffractive 
cross section of 80 mb, the LHC will produce on average 8 • 10^ inelastic p-p 
events per second, creating an extremely hostile radiation environment. For 
radiation studies, the bunch structure of LHC is not significant. The only 
scaling parameter for doses and fiuences is the inelastic interaction rate. 

The central tracker is expected to be exposed to a primary particle flux from 
the interaction region, and the main concern is radiation damage of silicon de- 
tectors. Without loss of generality, the radiation field simulated for the CMS 
silicon tracker geometry is considered in the following calculations [17]. The 
high magnetic field imposes a pr cut-off on charged particles, so that a signif- 
icant proportion of the most damaging low energy particles never reach the 
outer tracker layers. In addition, the average kinetic energy rises with increas- 
ing pseudorapidity. The spectra of charged hadrons (pions, kaons and protons) 
simulated for the positions of the silicon layers are taken from Reference [17]. 
The flux decreases by a factor 50 going from r = 20 cm to r = 100 cm, and 
most of the low energy particles disappear. The hadrons are predominantly 
low-energy charged pions and protons, that are present in different amounts 
and have different energy spectra as a function of the distance in respect to 
the interaction point, and of the pseudorapidity. In all cases, the pions are the 
dominant particles. 

In figure 1 the energetic differential generation rate of defects is presented. 
Each spectrum has been obtained as a convolution of the simulated hadron 
flux in the tracker cavity for CMS experiment and the energetic dependence 
of concentration of primary defects (CPD) for protons and pions in the same 
energy range. 

The calculations have been performed for two extreme positions in the tracker 
cavity: (1): r = 20 cm, z = ^ 140 cm, and (2): r = 100 cm, z = 140 H- 280 
cm. 

In the figure, the area under each curve represents the integral generation rate 
of vacancy-interstitials pairs. The values obtained for the first position (1) 
are: 6.2 • 10^ VI pairs/cm^/s for pions and 5.6 • 10'' VI pairs/cm^/s for protons, 
while for the second one (2) these are 8.1 • 10^ VI pairs/cm^/s and 3.1 • 10^ VI 
pairs/cm^ /s for pions and protons respectively. 
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Fig. 1. Differential energetic generation of vacancy-interstitial pairs for pions and 
protons, for two positions (1): r = 20 cm, z = 140 cm, and (2): r = 100 cm, 
z = 140 -\r 280 cm. in the tracking cavity of LHC. 



We would like to mention that the main contribution for protons comes from 
the lowest energy region, while for pions the maximum shifts from around 200 
MeV to 800 MeV passing from position (1) to position (2). For pions, the 
CPD dependence is cut at 20 MeV. 

The second type of application discussed in the present paper refers to the 
radiation field produced by cosmic rays. From these particles, the most im- 
portant contribution comes from protons. The primary proton spectrum in 
the kinetic range 0.2 to 200 GeV, in the neighbourhood of the Earth, at an 
altitude of about 380 km, was measured by the Alpha Magnetic Spectrometer 
(AMS) during space shuttle flight STS-91. The complete data set combining 
three shuttle altitudes and including all known systematic effects is given in 
reference [18]. The convolution of this spectrum with the CPD for protons 
is presented in Figure 2, and it corresponds to a generation rate of vacancy- 
interstitial pairs of 2 • 10^ VI pairs/cm^/s. This represents a nearly seven order 
of magnitude lower generation rate in respect to the higher one calculated for 
the tracking cavity of the LHC experiments. 
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Fig. 2. Differential energetic generation rate of vacancy-interstitial pairs by cosmic 
protons in the space near the Earth. 

4 Estimated damage in silicon detectors 



Silicon used in high energy physics detectors is n-type high resistivity (1-^6 
Kfl- cm) phosphorus doped FZ material. In the last decade a lot of studies 
have been performed to investigate the influence of different impurities, espe- 
cially oxygen and carbon, as possible ways to enhance the radiation hardness 
of silicon for detectors in the future generation of experiments in high en- 
ergy physics - see, e.g. references [19,20]. Some peoples consider that these 
impurities added to the silicon bulk modify the formation of electrically ac- 
tive defects, thus controlling the macroscopic device parameters. The effect of 
oxygen in irradiated silicon has been a subject of intensive studies in remote 
past. Empirically, it is considered that if the silicon is enriched in oxygen, 
the capture of radiation generated vacancies is favoured by the production of 
pseudo-acceptor complex vacancy-oxygen. Interstitial oxygen acts as a sink 
for vacancies, thus reducing the probability of formation of divacancy related 
complexes, associated with deeper levels inside the gap. These conclusions are 
confirmed by the present model. 

The concentrations of interstitial oxygen and substitutional carbon are silicon 
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are strongly dependent on the growth technique. In high purity Float Zone 
Si, oxygen interstitial concentrations are around lO'^^ cm~'^, while in the oxy- 
genation technique developed at BNL, an interstitial oxygen concentration 
of the order 5 • 10^^ cm~^ is obtained. These materials can be enriched in 
substitutional carbon up to [Ci\ 1.8 • 10^^ cm~'^. 



4-1 Changes in the microscopic material properties 

In Figure 3 a) to f), the formation and time evolution of the vacancy-oxygen, 
vacancy-phosphorus, divacancy, divacancy-oxygen, carbon interstitial-oxygen 
interstitial and carbon interstitial-carbon substitutional are presented. Two 
types of silicon have been considered: the "standard" material, containing the 
following impurities concentrations: 10^^ cm~^ atoms of phosphorus, 2 • 10^^ 
cm~^ atoms of oxygen, and 5 • 10^^ cm~'^ atoms of carbon; and the "oxygened" 
one, containing 10^^ cm~^ atoms of phosphorus, 4 ■ 10^^ cm""^ atoms of oxy- 
gen, and 5 • 10^^ cm~^ atoms of carbon respectively, and for the two generation 
rates of VI - pairs discussed above, namely 7 ■ 10^ and 2 • 10^ VI pairs/cm^/s. 
The curves in the figures are labelled as follows: (1): standard silicon, LHC 
irradiation rate, (2): oxygenated silicon, LHC irradiation; (3): standard sili- 
con, cosmic irradiation rate; and (4): oxygenated silicon, cosmic exposure. As 
could be observed from the figure, the content of oxygen in silicon influences 
especially defects formation in the case of high rates of generation of vacancy- 
interstitial pairs. The increase of the initial oxygen concentration in silicon, 
conduces, after ten years of operation in the LHC environment, characterised 
by a high and constant generation rate, to the increase of the concentrations 
of VO and CjOj centres, and to the decrease of the concentrations of V2, VP 
and CiCs ones. With the increase of oxygen concentration, an increase of the 
V2O generation rate is observed. It is interesting to observe that in almost 
all cases, an equilibrium in reached between generation and annealing, and a 
plateau is obtained in the time dependence of the concentrations. The slowest 
is, in this respect, V2O, that has the highest binding energy. 

As underlined above, vacancy-oxygen formation in oxygen, enriched silicon is 
favoured in respect to the generation of V2, V2O and VP centres. At high oxy- 
gen concentrations, the concentrations of VO centres attain a plateau during 
the 10 years period considered. 

After cosmic proton irradiation, the effects are strongly different. For this gen- 
eration rate, the increase of the oxygen concentration produces the decrease of 
the concentration of all centres, with the exception of the VO concentration, 
that, at these rates, is not influenced by the oxygen content, and of the CiOi 
concentration, where an increase is observed. As a consequence of the small 
rate of generation of vacancy-interstitial pairs, after ten years of operation. 
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Fig. 3. Time dependence of the concentrations of a): VO, b): VP, c): V2, d): V2O, 
e): CiOi and f): CiCg induced in standard and oxygenated silicon, in conditions of 
LHC and cosmic continuous irradiation rates (see text). 



the equilibrium between generation and annealing is not reached, the concen- 
trations of defects being, with the exception of VP (that has a relatively low 
binding energy), slightly increasing functions of time. 

All the processes have been calculated for 20^0 temperature. Thermal gener- 
ation has been taken into account in all cases, although it is important only 
for the silicon exposed to cosmic protons. 
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Also, it is necessary to emphasise the importance of the irradiation and an- 
neahng conditions (initial material parameters, type of irradiation particles, 
energetic incident particle spectra and their flux, temperature) on defect evo- 
lution. These aspects have been discussed in previous papers [13,16]. 

4.2 Macroscopic modifications of detector parameters - the leakage current 

The dark current of a reverse biased p — n junction is composed of the fol- 
lowing terms: the drift current, due to the drift of minority carriers, the gen- 
eration current, due to carrier generation on the midgap energy levels inside 
the depleted region and surface and perimetral currents, dependent on the 
environmental conditions of the surface and the perimeter of the diode. The 
formation, during and after irradiation, of defects with associated energy lev- 
els inside the gap conduces to the increase of the generation current, since the 
ease with which a mobile carrier can traverse the gap is greatly enhanced by 
intermediate levels. 

Inside the depleted zone, n,p <^nt {ut is the intrinsic free carrier concentra- 
tion); each defect with a bulk concentration Nt causes a generation current 
per unit volume of the form [21]: 

I -qU -q<vt> n,^^^^^(^^_^^)/^^^ ^ [14) 

where 7„ and 7^ are degeneration factors, (T„ (up) are the cross sections for 
majority (minority) carriers of the trap, Ei = [Ec — Ey)/2 and < > is the 
average between electron and hole thermal velocities. In the S ho ckley- Read- 
Hall model used for the calculation of the reverse current, each defect has one 
level in the gap, and the defect levels are uncoupled, thus the current is simply 
the sum of the contributions of different defects. Two parameters characteris- 
ing the defects enter in the calculation of the generation current: their energy 
position in respect to the intrinsic level, and their cross section. Only near 
midgap energy levels are important, and in this paper the contributions com- 
ing from V2O, V2 and VP have been taken into account. An average between 
the values of the energy levels and cross sections reported in the literature (see 
compilations [5,14]) for V2 and VP have been introduced in the calculations, 
and averaged while for V2O this is true only for the energy level. In the lack of 
reported data for the cross section, for the V2O centre, the value 10"^^ cm~^ 
has been used. 

Divacancy has three energy levels in the band gap. The equilibrium statistics 
for this case is formally different for this case from that for the same number 
of independent levels, since the occupancy of the levels is now interdependent. 
In the present calculations, the interaction of the V2 energy levels has been 
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neglected (in accord with Sah and Schotkley affirmation [22] that the inde- 
pendent and interacting energy level cases are indistinguishable if the energy 
levels are more than a few ksT apart). 

In Figure 4, the separate contributions of V2O, V2 and VP defects to leakage 
current for "standard" and "oxygenated" silicon, in the irradiation conditions 
supposed for LHC and in the cosmic near Earth orbits, are represented. The 
calculations have been performed for 20°C, under continuous irradiation. The 
bands representing the maximal uncertainties due to the energy position of 
the defects and to their cross sections are also drawn as dotted lines. On can 
observe that in the conditions of LHC generation rate, the values of the total 
current are nearly the same for standard and oxygenated silicon: the higher 
contributions from the V2 and VP centres (standard silicon) are counterbal- 
anced by the increase of the V2O concentration (oxygenated silicon), that is 
the nearest to the midgap. For smaller generation rates, oxygenated silicon is 
a better choice from the point of view of the leakage current, as could be seen 
from the case of cosmic irradiation. 



5 Summary 

The phenomenological model developed previously to explain defect genera- 
tion and evolution in silicon has been used to evaluate the damage induced 
by hadron fields, for two classes of applications in high energy physics: at the 
new generation of colliders and in space applications. The generation rates of 
vacancy-interstitial pairs have been calculated as convolutions of the hadron 
spectra with the energy dependencies of the CPD. The time dependence of 
the concentrations of stable defects has been calculated for "standard" and 
"oxygenated" silicon, in conditions of continuous irradiation during 10 years, 
at 20°C, for generation rates of vacancy interstitial pairs corresponding to the 
applications mentioned before. The increase of the oxygen content in irradi- 
ated silicon conduces to the increase of the concentrations of oxygen related 
defects, VO, CiOi and V2O, and to the decrease of the V2, VP and CiCg ones. 
For long term operation and high generation rates of VI pairs, comparable 
leakage currents are expected in standard and oxygenated silicon p-n junc- 
tions. The beneficial influence of a higher oxygen content in silicon becomes 
visible with the decrease of the generation rate. 
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